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-»^ ■ Abstract 

43, 
Qh! We predict a new mechanism for the spin light of neutrino (SLv) that can be 

emitted by a neutrino moving in gravitational fields. This effect is studied on the 

basis of the quasiclassical equation for the neutrino spin evolution in a gravitational 

field. It is shown that the gravitational field of a rotating object, in the weak-field 

limit, can be considered as an axial vector external field which induces the neutrino 

^ spin procession. The corresponding probability of the neutrino spin oscillations in the 

gravitational field has been derived for the first time. The considered in this paper 

SLv can be produced in the neutrino spin-flip transitions in gravitational fields. It is 

shown that the total power of this radiation is proportional to the neutrino gamma 

factor to the fourth power, and the emitted photon energy, for the case of an ultra 

relativistic neutrino, could span up to gamma-rays. We investigate the SLv caused by 

both gravitational and electromagnetic fields, also accounting for effects of arbitrary 

moving and polarized matter, in various astrophysical environments. In particular, 

we discuss the SLv emitted by a neutrino moving in the vicinity of a rotating neutron 

star, black hole surrounded by dense matter, as well as by a neutrino propagating in 

the relativistic jet from a quasar. 

1 Introduction 

It is well known that a neutrino with non-zero mass has non-trivial electromagnetic 
properties [1] (for the recent study of the electromagnetic properties of a massive neutrino 



see [2,3]). In particular, the Dirac massive neutrino possesses non- vanishing magnetic 
moment. Even the massive Majorana neutrino can have transitional magnetic moment. 

In a series of our papers we have developed the Lorentz invariant approach to neutrino 
spin (and also flavour [4,5]) evolution in different environments accounting for the presence 
of not only electromagnetic fields [6,7] and weakly interacting with neutrino matter [8] but 
also for other types of neutrino non- derivative interactions [9]. On this basis we have 
also predicted [10] the possibility for a new mechanism of electromagnetic radiation by 
neutrino moving in background matter and/or electromagnetic fields. The new mechanism 
of electromagnetic radiation ( we have named [10] this radiation "spin light of neutrino" 
(SLv) ) originates from the neutrino spin precession that can be produced whether by weak 
interactions with matter or by electromagnetic interactions with external electromagnetic 
fields. The latter possibility was also considered before in [11]. A review and some new 
results on our studies of neutrino oscillations in matter and external fields is given in [12]. 

As it is shown in [10], the total power of the SLv is not washed out even when the 
emitted photon refractive index in the background matter is equal to unit. That is why 
the SLv can not be considered as the neutrino Cherenkov radiation (see, for example, [13] 
and references therein). The total power of the SLv in matter is proportional to the 
fourth power of the matter density and the neutrino Lorentz factor. The SLv is strongly 
beamed in the direction of neutrino propagation and is confined within a small cone given 
by d ~ 7 _1 . The energy range of emitted SLv photons, for the case of the relativistic 
neutrino, could span up to gamma-rays. These properties of SLv enables us to predict 
that this radiation should be important in different astrophysical environments (quasars, 
gamma-ray bursts etc) and in dense plasma of the early Universe. 

In this paper we should like to introduce a new type of the SLv that could originate 
from the neutrino spin precession induced by a gravitational field in different astrophys- 
ical and cosmological settings. This study is based on consideration of the neutrino spin 
evolution problem in gravitational fields. We also consider the neutrino spin oscillations in 
gravitational fields and derive for the first time the corresponding probability. 

The problem of the neutrino flavour oscillations and neutrino spin evolution in the 
presence of gravitational fields have been studied in many papers. The gravitational effect 
on flavour neutrinos from collapsing stars was considered in [14] . Neutrino spin procession 
and oscillations in gravitational fields were also studied in [15-20]. The important obtained 
result of these studies is that if the contribution of gravitational interaction is diagonal in 
spin space (that is the case of gravitational fields of non-rotating massive objects), then 
gravity cannot produce the neutrino spin procession on its own. As it is pointed out 
in [16], the off-diagonal terms in neutrino spin space could appear from the interaction of 
a neutrino magnetic moment with a magnetic field. 

A rather detailed discussion on various aspects of the neutrino spin and chiral dynamics 
in the presence of gravitational fields can be found in the recent papers [21-25] where, in 
particular, the case of a strong field in a Schwarzschild space-time background have been 
also considered. The case of a weak field was also studied in [18] . However, the possibility 
of the electromagnetic radiation due to the neutrino spin procession in a gravitational field 
has never been discussed before. 



Here below we apply the Lorentz invariant approach, based on the generalization of the 
quasiclassical Bargmann-Michel-Telegdi equation for the case of different non-derivative 
interactions of neutrino with external fields [6,8,9], to the study of the neutrino spin 
evolution in weak gravitational fields. We found out that weak gravitational fields of a 
rotating source enters the neutrino spin evolution equation as an axial- vector field, that, 
as it has been shown in [9], produces the spin precession. Then, using the main idea 
of [10] that a neutrino with spin processing in the external environment should radiate 
electromagnetic waves, we investigate this spin light of neutrino in gravitational fields in 
the vicinity of the rotating neutron star, the accretion disk of a rotating black hole and in 
the relativistic jet of a quasar. 

2 Neutrino spin evolution equation in gravitational 
fields 

Following the study of ref. [18], our starting point is the Dirac equation for a neutrino 
in a curved space-time [26]: 

[i^(x)D l , - m]* = 0, (1) 

where D^ = d^ + T^ is the covariant derivative operator, T^ = — {i/S)e u a e v i ) ;n\l a 1 l b } is the 
spin connection, and e^ are the tetrads. Latin indices refer to a local Minkowski frame, 
while greek indices refer to general curvilinear coordinates. The gamma matrices 7 M = 7 a e^ 
satisfy the conditions {7^, 7^} = 2g^ u and D^j u = 0. The metric tensor g^ v and the tetrads 
e^ are related by 

9^ = Vabefcl, (2) 

where r) a b is the Minkowski metric. 

We consider below the limit of a weak gravitational field which is appropriate for the 
most of astrophysical systems. Indeed, we are also restricted to the weak-field limit by 
our Lorentz invariant approach in description of the neutrino spin evolution on the basis 
of the generalized Bargmann-Michel-Telegdi equation. Therefore, we perform our analysis 
within the linearized theory of gravity by assuming that the metric can be written as 

and treat the second term in the right-hand side as a perturbation over the flat Minkowski 
space. In this approach h^ u is the gravitation field, and the quantity k is related to the 

Newton's constant Gn, 

k = ^8irG N . (4) 

As it has been shown in [18], the corresponding Hamiltonian for the considered linear 
approximation is 

H = 7 °m(l - k h 00 ) + op- ^{h 00 , ap} + - \{ti\ a i p>} + + k {h, p} + + k ([V x h] s) , (5) 



where p = — iftV, ex = 7°7, and s = |S = |7°7 5 7- Note that in the limit of a weak 
gravitational field it is possible to treat gravitational perturbations h m = h as an external 
field contributing to the neutrino interaction Lagrangian. Starting from the neutrino inter- 
action Lagrangian C accounting for general types of neutrino non-derivative interactions 
with external fields, 

-C =g s s(x)uu + g p 7i-(x)vj 5 v + g^^^v^^v + g a N l {x)v^ iX 'fv 

+ Ir^Ptvi/ + |n^pa^75^, (6) 

where s, ir, V^ = (V f) ,V),A fi = (A , A), T^ = (a, b),II MV = (c, d) are the scalar, pseu- 
doscalar, vector, axial-vector, tensor, pseudotensor fields, respectively, and Qi (i = a, p, 
v, a, t, t') are the coupling constants, a^ v = (i/2)(j^y u — 7^7^), we get [9] the following 
expression for the corresponding Hamiltonian: 

U =9ssp 3 - ig p 7Tp 2 + g v (V° - (aV)) - g a (piA° - (SA)) 

- g t (p 3 CEb) + p 2 (Sa)) - ig' t {p^c) - p 2 (Sd)), (7) 

where a = 7°7, p x = -75, p 3 = 7 , p 2 = ipip?,. 
The last term in equation (5), 

ftspin = k( [V x h]s), (8) 

which is responsible for the neutrino spin interaction with a gravitational field, can be 
identified with an axial-vector part of the general Hamiltonian given by (7) if we define 
the axial- vector field A^ = (0, A) by the following relation 

<? a A = ^[Vxh], (9) 

where g a is the coupling constant of this field to a neutrino. Following the performed 
in [9] derivation of the quasiclassical spin evolution equation accounting now only for the 
axial- vector interaction given by (9), we get the evolution equation for the neutrino spin 
vector C, u in the considered gravitational field 

f^[CxG], (10) 

where t is time in the laboratory frame, and, as it follows from eq.(8) of [9], the vector G 

is 

k f r __ 7 



G = -- |[V x h] + T _^/3(/3[V x h])| . (11) 

Here /3 is the speed , 7 = E u /m u is the Lorentz factor of the neutrino. 

It should be noted that the derived neutrino spin evolution equation (10) accounts only 
for the gravitational field, however the possible effects of the neutrino interactions with 



the background matter and electromagnetic fields are not included here. The analogous 
neutrino spin evolution equation in the presence of electromagnetic fields and matter in 
the Minkowski space ( see refs. ( [6,9,10,12]) is 



B f 



dt 

7(B_l + -B|i- 

7 



M 0|l = 7/3 
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(13) 



(14) 



(15) 



(16) 



where F^ and Fn (F = B, E) are transversal and longitudinal (with respect to the direction 
of neutrino motion n = (3/(3) electromagnetic field components in the laboratory frame. 
For simplicity we neglect here the neutrino electric dipole moment, e = 0, and also consider 
the case when matter is composed of only one type of fermions (electrons). The general 
case of e 7^ and matter composed of different types of leptons is discussed in our papers 
mentioned above ( [6,9,10,12]). Note that n 



n e y 1 — fg in (15) and (16) is the invariant 
number density of matter given in the reference frame for which the total speed of matter is 
zero. The vectors v e , and £ e (0 ^ \( e \ 2 ^ 1) denote, respectively, the speed of the reference 
frame in which the mean momentum of matter (electrons) is zero, and the mean value 
of the polarization vector of the background electrons in the above mentioned reference 
frame. The coefficients p e ' are calculated if the neutrino interaction with matter is given. 
Therefore, within the extended standard model supplied with S77(2)-singlet right-handed 
neutrino z/r, we have 



P, 



(i) 



G, 



2V2fi 



P 



(2) 



G, 



2y/2fi 



(17) 



where Gf = Gf(1 + 4sin 2 6*jy). 

From the neutrino evolution equation (12) we have predicted that the neutrino spin 
precession can be induced not only by external electromagnetic fields but also by the 
neutrino weak interaction with particles of the background matter. Now from the neutrino 
spin evolution equation (10) in the gravitational field it follows that the neutrino spin 
procession can appear if the vector G is at least not zero. Thus, if the spin is processing 
the off-diagonal metric components h° l have to depend on the space coordinates. The 



appropriate example is provided by the Kerr geometry that corresponds to the gravitational 
field of a rotating with angular momentum L object. In this case 

h = A [Lxr]! ( 1 8) 

and the axial- vector field A can be written as 

ft ,A = G„ 3r(L ^-' JL . (19) 

Thus, we conclude that, as it follows from the neutrino spin evolution equation (10), 
the spin procession can appear if the neutrino is moving in a weak gravitational field of a 
rotating object. 

3 Probability of neutrino spin oscillations in gravita- 
tional fields 

Using the results of the previous section, it is possible to get the corresponding neu- 
trino spin oscillations Hamiltonian which can be used in derivation of the neutrino spin 
oscillations probability in the presence of the gravitational field. 

Let us consider the case when the neutrino is moving along the radial direction and 
decompose the vector G following 



G = G||+G_l, (20) 

where 

r-3 



G|| = -7^^nc()s^. (21) 



and 

G ± = ^n ± sin0 (22) 

2 r A 

are the longitudinal and transversal (with respect to the neutrino momentum) components 

of the vector G, n^ is the unite orthogonal to n and laying in one plane with n and L 

vector. In order to get the probability of the neutrino spin oscillations in the gravitational 

field we introduce the unit vector k 

"-§■ 

and rewrite eq.(10) in the form 

dC, v a 

dt r 

where with the use of (20), (21), and (22) we have 

a = \ A'y 2 cos 2 9 + sin 2 9. 

7 v 



3 [Cxk], (23) 



If we expand the neutrino spin vector £„ over the basis determined by the vector G, so 
that 



where (C^k) = and (GCt) — 0, then the solution of eq.(23) in terms of £jl, Cx and Cy 
(C* = y/Cx + Q), is given by 



Cx = -\/i-Co sin 



O 



C = Co- 



2(5r 2 
a 



2[3r 2 J ' 



(24) 



Here Co is a constant determined by the initial conditions. Finally, for the probability of 
the neutrino oscillations in the gravitational field, described by the vector G, we get 



Pu L ^u R (t) = ^[l + (Cn) 



(25) 



were C, v is the solution of eq.(23) given by (24). 



4 Spin light of neutrino in gravitational fields 

From our previous studies [10] of the SLu we know that this radiation is emitted if the 
neutrino magnetic moment is processing, no matter what is the cause of this procession. 
Therefore, using the results for the SLu in matter and electromagnetic fields and the 
analogy between neutrino spin evolution eqs. (10) and (12) we get for the total power of 
the spin light of neutrino in the gravitational field 



T -1? 2 



4(G 2 ) 2 + G 2 



(26) 



where \x is the neutrino magnetic moment. 

Let us now consider the new phenomenon of the SLu in the case of neutrinos moving in 
the gravitational fields of a rotating neutron star along its radial direction. From eqs. (20), 
(21), and (22) it is easy to see that for the ultra relativistic neutrinos ( for the neutrinos 
from neutron stars the Lorentz factor could be of the order of 7 ~ 10 9 ) the total power of 
the SLu is maximal when neutrinos move along or opposite to the angular momentum L: 



t gT 



(~<2 r2 



1-^ 



r<2 7-2 

9r 4 



(27) 



For a neutron star with the radius ~ 10 km we get 



2_,4 



I gr = 3/i7 



1 + 



1 ( \ km x 



144 V r 



^ (28) 



It follows that the total radiation power of the SLv in the gravitational field is proportional 
to ~ 7 4 and the main contribution is given by the second term in (27) which originates 
from the derivative term in the general expression (26). It is also possible to show that 
there is a strong beaming effect and that the radiation is confined within a small cone in 
the direction of the neutrino propagation given by the angle 59 ~ 7" 1 . Thus, we predict 
that the angular distribution of the SLv in the gravitational field of a rotating neutron 
star is not isotropic even if neutrinos are moving symmetrically in all radial directions from 
the neutron star. In this case the SLv in the gravitational field is radiated more effectively 
by neutrinos moving along or against the axis of the neutron star rotation. 

It is also possible to consider the combining effect of the electromagnetic, weak and 
gravitational fields on neutrino spin procession. From (10) and (12) we get for the neutrino 
spin evolution in this case 

dC 2 

■f = -[Cx (//Bo + /iMo + G)] , (29) 

at 7 

where Bo, Mo are given by eqs.(13) and (14) respectively. The gravity effect enters through 
the vector G, which in the considered case of a rotating neutron star, is given by eq.(20). 

It is interesting to compare the contributions to the SLv from the neutrino interaction 
with the background matter and with the background gravitational field. For illustration, 
let us investigate the two cases: i) the SLv from a rotating black hole produced by neutrinos 
moving in accretion disk, and ii) the SLv from a quasar produced by neutrinos moving 
along the relativistic jet. We take matter consisting of hydrogen. In the first case we assume 
that matter in the accretion disk is of constant density and slowly moving along the radial 
direction. In the second case we consider matter of the jet moving with relativistic speed. 
We omit here the electromagnetic interaction but account for the neutrino couplings with 
matter and gravitational field, aiming to consider also the interplay between them. This 
modelling of the background environments, although is rather simplified, could reveal the 
main features of the effects under consideration. 

We start with the general expression [9] for the neutrino spin evolution Hamiltonian 
and extract for the further consideration the terms accounting for the neutrino interactions 
with the gravitational field and matter: 

U = 9a (EA) + ( 75 /°(x) + (S f (*))) , (30) 

where g a A is given by (9). The matter term contribution in (30) is determined in the 
general form as 

f=e,p,n 



qf } = (/<{» -2Q"> sin 2 6 W + S ef ), qf = -(/<» +S„), (32) 

[0, / = n, p, 

however in the considered case of matter composed of hydrogen there is no contributions 
from neutrons n. In this equations, I^ L denotes the third component of weak isospin for 
the fermion /, and Q^ is the electric charge of the fermion. The matter current and 
polarization are [8], respectively, 

f = (n f ,n f v f ), (33) 

-/C/v / ,nK /v / r^+^4= ) J=e > p - ^ 

l + \J l - v f) 
Finally, for the neutrino spin evolution equation in the considered case we have 

^ = - [C x ( M M + G)] , (35) 

at 7 

where the matter term according to eqs.(30), (31), and (32) is 

M = 7/3(V- T ^(W)-f, (36) 

and the gravitational field effect is described by the vector G which is given by (20). Thus, 
for the total radiation power of the SLv in matter and the gravitational filed we find, 

/gr.matt = "gV ^((//Mq + G) 2 ) 2 + (/lM + G) 2 J . (37) 

It is evident that the total radiation power, I grima tt, is composed of the three contributions, 

-'gr,matt -'gr "T~ ^matt ~r -'gr+matt; V^"J 

where I gT is the radiation power due to the neutrino magnetic moment interaction with the 
gravitational field, I ma tt is the radiation power due to the neutrino weak interaction with 
the background matter, and I gr + ma tt stands for the interference effect of the gravitational 
and weak interactions. In an analogy with the neutrino spin flip in matter (see [10]) we 
predict the similar effect for neutrino spin in gravitational fields: the initially unpolarized 
neutrino beam (equal mixture of active left-handed and sterile right-handed neutrinos) can 
be converted to the totally polarized beam composed of only vr. 

Let us turn to the discussion on the SLv in the two particular cases mentioned above: 

i) It is obvious that in this case at small distances r from the center of the black hole 

the contribution to the SLv coming from the gravitational interaction dominates over 

9 




Figure 1: The total radiation power of the spin light of neutrino, I(r) = I gT + I m att + 
-fgr+matt, versus radial distance r is shown by the red solid line. The blue dashed and green 
dash-dotted lines correspond, respectively, to the gravitational, I gr , and matter, / ma tt , 
contributions without account for the interference term / gr + m att- 

that from the interaction with matter. As neutrino propagates away from the center, 
the gravity contribution fades and starting from a certain distance from the black hole 
center the matter interaction contribution to the SLu becomes predominant. The total 
radiation power, / gr ,matt = I( r )i ver sus the distance r is shown in Fig. 1. The green dash- 
dotted line corresponds to the matter contribution, while the blue dashed line stands for 
the gravitational field contribution to the SLv. For the neutrino gamma factor, neutrino 
magnetic moment and matter density we take, respectively, the following values of 7 = 10 12 , 
\i = lCT 10 /io and n = 10 24 cm~ 3 , where /i is the Bohr magneton. It is supposed that the 
gravitational field is produced by the rotating object with the solar mass M = M , and 
the angular momentum is chosen to be equal to the maximal allowed value L = ^/(AG^) 
( see, for instance, [27]) where r is the Schwarzschild radius, 
ii) In this case, following eq.(35) the total radiation power is 



T _^1 2 4 

'gr,matt „ p 7 



GatL 



Gp ^-v 

f=e,P 



n 



A v / 



(3) 



9 G 2 N f3 2 L 2 

4 7"8 



(39) 



If matter is moving with relativistic speed along the direction of the neutrino propagation 
then the matter contribution to / gr ,matt is washed out because of the presence of the term 
(v/ — (3) ~ . This happens due to the particular properties of the neutrino oscillations 



10 



in moving matter (see [4,5,8,9,12]). Thus, we conclude that in the case of the neutrinos 
moving in the relativistic jet of a quasar the SLv could be produced by the gravitational 
field, whereas there is no important contribution from interaction with matter. 

5 SLv photon energy in gravitational fields 

Now we discuss the energy range of the SLv photons emitted in a background where 
only the presence of the gravitational field of a rotating object generates the neutrino spin 
procession (it is supposed that neither matter nor electromagnetic fields of the background 
gives important contribution to the neutrino spin evolution). In order to get an estimation 
for the characteristic scale of the energy of the emitted photons we suppose that the 
variation of the vector G (which describes the influence of the gravitational field) on 
the neutrino travelled distance Ar is much less then the frequency of the neutrino spin 
precession, 

AG 

-— < uAr. (40) 

If we consider the neutrino propagating along the rotation axis, 9 = 0, and also take into 
account followed from eqs.(21) and (22) relation G\\ ^> G± then the emitted photon energy 
in the rest frame of the neutrino is 

To calculate the photon energy in the laboratory reference frame we should use the well 
known formula 

<" = ^o^P\, (42) 

1 — p COS V 

where d is the angle between the neutrino speed and the direction towards the observer. 

Due to the fact that the SLv is strongly beamed within a small cone there is a reason to 

consider the case of d ~ 0. Therefore, the emitted photon energy in the laboratory frame 

is 

GnL 2 

uj ~ cj 7 ~ — —7 . (43 J 

For the set of parameters that we have just considered in the case of a rotating black hole, 
eq.(43) reads 

uo ~ KT 11 x 7 2 (—) eV. (44) 

For 7 ~ 10 12 and r ~ 10 tq we obtain uo ~ 10 GeV. Note that for the mass of neutrino 
m u ~ 1 eV we also obtain (u/E u ) ~ 10~ 2 ^C 1, i.e. the quasiclassical approach to the 
neutrino spin evolution is valid in this case. These properties of SLv enable us to predict 
that this radiation should be important in different astrophysical environments (quasars, 
gamma-ray bursts etc) and in dense plasma of the early Universe. 
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6 Conclusion 

We also should like to point out that although the predicted above and in [10] the 
spin light of neutrino in matter, electromagnetic and gravitational fields is supposed to be 
radiated in the process 

v\ -»• v 2 + 7, (45) 

without change of the neutrino flavour state (the neutrinos in the initial and final states, 
v\ and 1S2, are of the same flavour), it is possible to generalize, with a minor modification, 
the corresponding equations for the case when the SLv is radiated due to the neutrino 
transition magnetic moment interaction with the background fields. The latter corresponds 
to models of the Dirac neutrinos with non-diagonal magnetic moments or the Majorana 
neutrinos which could also have transitional (magnetic) moments. 

In conclusion, we consider the neutrino spin evolution in presence of gravitational fields 
and derive for the first time the corresponding neutrino oscillations probability. We also 
predict the possibility of a new mechanism of the spin light of neutrino [SLv) in a gravita- 
tional field which, as we expect, could have consequences in astrophysical and cosmological 
settings. As examples, we analyzed the SLv in particular cases of neutrino moving in the 
vicinity of a rotating neutron star, black hole and the relativistic jet of a quasar. 
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